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ABSTRACT 

The stellar population in the Galactic halo is characterised by a large fraction of carbon-enhanced metal-poor (CEMP) stars. Most 
CEMP stars have enhanced abundances of i-process elements (CEMP-i stars), and some of these are also enriched in r-process 
elements (CEMP-i/r stars). One formation scenario proposed for CEMP stars invokes wind mass transfer in the past from a thermally- 
pulsing asymptotic giant branch (AGB) primary star to a less massive companion star which is presently observed. In this work we 
generate synthetic populations of binary stars at metallicity Z = 0.0001 ([Fe/H] -2.3), with the aim of reproducing the observed 

fraction of CEMP stars in the halo. In addition, we aim to constrain our model of the wind mass-transfer process, in particular 
the wind-accretion efficiency and angular-momentum loss, and investigate under which conditions our model populations reproduce 
observed distributions of element abundances. We compare the CEMP fractions determined from our synthetic populations and the 
abundance distributions of many elements with observations. Several physical parameters of the binary stellar population of the halo 
are uncertain, in particular the initial mass function, the mass-ratio distribution, the orbital-period distribution, and the binary fraction. 
We vary the assumptions in our model about these parameters, as well as the wind mass-transfer process, and study the consequent 
variations of our synthetic CEMP population. The CEMP fractions calculated in our synthetic populations vary between 7% and 
17%, a range consistent with the CEMP fractions among very metal-poor stars recently derived from the SDSS/SEGUE data sample. 
The resulting fractions are more than a factor of three higher than those determined with default assumptions in previous population- 
synthesis studies, which typically underestimated the observed CEMP fraction. We find that most CEMP stars in our simulations 
are formed in binary systems with periods longer than 10,000 days. Pew CEMP stars have measured orbital periods, but all that do 
have periods up to a few thousand days. Our results are consistent only if this small subpopulation represents the short-period tail 
of the underlying period distribution. The results of our comparison between the modelled and observed abundance distributions are 
significantly different for CEMP-j/r stars and for CEMP-i stars that are not strongly enriched in r-process elements. Por these stars, 
our simulations qualitatively reproduce the observed distributions of carbon, sodium, and heavy elements such as strontium, barium, 
europium, and lead. Contrarily, for CEMP-^/r stars our model cannot reproduce the large abundances of neutron-rich elements such as 
barium, europium, and lead. This result is consistent with previous studies, and suggests that CEMP-i/r stars experienced a different 
nucleosynthesis history to CEMP-i stars. 
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1. Introduction 

The population of old metal-poor stars observed in the Galactic 
halo carries information about the conditions under which the 


Milky Way was formed. In the past three decades, the HK sur- 

vey ([Beers et al. 

1985[[1992[|, the Hamburg/ESO survey (HES; 

[Christlieb et al. 

2001 200811, and the SDSS/SEGUE survey 

(York et al. 2000 

Gunn et al.|2006[ Yanny et al.|2009 1 have col- 


nomenclature proposed by Beers & Christlieb ( 2005 | l, CEMP 
stars have [C/Ee] > H-1. More recently a number of authors have 
adopted the criterion [C/Ee] > -1-0.7 to define CEMP stars (e.g. 
Aoki et al.|2007 Lee et al. 2013[[^ng et al. 2013|l. The fraction 


of CEMP stars among metal-poor stars is observed to rise with 
increasing distance from the Galactic plane (Erebel et al.|2006 


Carollo et al. 2012|l. A strong increase in the cumulative fraction 


lected spectra of a large sample of v^ metal-poor stars (VMP, 
that is, stars with abundances of irorQ[Pe/H] < -2 ), and have 
made it possible to investigate their dynamical and chemical 
properties. These studies reveal a large proportion of carbon- 
enhanced metal-poor (CEMP) stars. According to the original 


of CEMP stars with decreasing metallicity has been confirmed 
by many studies, from about 3% at [Ee/H] < -1 to around 75% 
at [Ee/H] < -4 (e.g. [Cohen et al.|2005| [Marsteller et al.|2005t 
Erebel et al. 20061|Lucatello et al.|20061|Carollo et al.|2012l|Aoki| 

et al.|201^ ^ee et al.|2013[ Yong et al.|2013 Placco et al. 2014| l. 


' The relative abundance of two elements X and Y is [X/Y] = 
logjofAx/Ay) - logiofAx/Alylo, where Nx and Ay refer to the num¬ 
ber density of atoms of X and Y, respectively, and G indicates the solar 
values. 


CEMP stars are usually classified into four categories ac¬ 
cording to the observed abundances of barium and europium, 
two elements associated with the slow (s-) and rapid (r-) 
neutron-capture process, respectively. 
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CEMP-s stars exhibit enhanced abundances of i-elements, 
and are defined by [Ba/Fe] > +0.5 and [Ba/Eu] > 0. CEMP- 
s stars account for at least 80% of all observed CEMP stars 
( |Aold et al.|2007| l. A number of authors adopt slightly differ¬ 
ent criteria, e.g. [Ba/Ee] > +1 and/or [Ba/Eu] > +0.5 (|Beers 


|& Christlieb|20(^[|Jonsell et al.|2006[|Masseron et al.|2010| l. 

• CEMP-s I r ifar.!|jare CEMP-i stars that also enriched in eu¬ 
ropium and r-elements, [Eu/Ee] > +1. 

• CEMP-r stars are a rare CEMP subclass enriched in ele¬ 
ments produced purely by the r-process, and are defined by 
the criteria [Eu/Fe] > +1 and [Ba/Eu] < 0. 

• CEMP-no stars do not exhibit peculiar enhancements in the 
elements heavier than iron ( jAoki et al.|20021 l. 

The formation scenario of CEMP stars is still uncertain, and 
different categories of CEMP stars are likely to have different 
formation channels. Several mechanisms have been proposed 
to explain the large carbon abundances, for example: (a) the 
gas cloud in which these stars were formed was already en¬ 
riched in carbon produced by zero-metallicity and/or rapidly 
rotating stars ([Mackey et ah 2003[ |Meynet et al. 2006[ |2010[ 
Maeder et al.|2014|l or exp^ed by the faint supernovae associ¬ 


ated with the first generation of stars (jUmeda & Nomoto|2003 


|2005[[Nomoto et al.|2013| ); and (b) the carbon-enhanced materia 
was accreted from the envelope of a thermally-pulsing asymp¬ 
totic giant branch (AGB) primary star in a binary system onto 
the presently observed low-mass companion. The mass-transfer 
scenario provides a natural explanation for the chemical abun¬ 
dances of CEMP-i stars because carbon and i-elements are pro¬ 
duced during AGB nucleosynthesis, and the observed fraction of 
CEMP-s stars with a binary companion is consistent with the hy¬ 
pothesis of all CEMP-s stars being in binary systems ( |Lucatello| 
et al.|200^ Starkenburg et al.|2014| l. By contrast, many CEMP- 


no stars show no indication of a binary companion ( Norris et al. 
|2013[ Hansen et al., in prep.), and they exhibit low abundances 
of s-elements, whereas large enhancements of s-elements are 
normally produced in AGB stellar models, except in some mod¬ 
els of mass M* > 3 Mq. The origin of CEMP-s/r stars is also an 
open issue. Current stellar models indicate that in AGB stars the 
density of free neutrons is not sufficiently large for the r-process 
to take place, hence the enhancement of r-elements in CEMP- 
s/r stars remains unexplained. 

A number of authors have simulated populations of very 
metal-poor stars with the aim of reproducing the observed frac¬ 
tion of CEMP stars. The largest uncertainties in these models 
include the amount of carbon and s-elements produced by AGB 
nucleosynthesis as a function of the stellar mass and metallicity, 
the shape of the initial mass function (IMF) of the early Galaxy, 
the binary fraction, the efficiency of the mass-transfer process, 
and the range of separations in which it occurs. [Luc atello et al. 


([2005 |l and Komiya et al. ( 2007|) argue that the IMF in the early 


Galaxy was weighted towards intermediate-mass stars to ac¬ 
count for the large fraction of CEMP stars at low metallicity. 
Suda et al. (2013|l and Lee et al. ([2014 1 suggest that the transi¬ 


tion between the early (top-heavy) and the present-day IMF oc¬ 
curred over the metallicity range between [Fe/H] = -2.5 and 
[Fe/H] = -1.5. However, [Pols et ah] P012l l show that, with 
an IMF biased towards intermediate-mass stars, many nitrogen- 
enhanced metal-poor stars are produced, whereas very few are 
observed. 

jlzzard et al.j ( |2009| l and [Abate et al. ( |20131 l model the 
CEMP population at metallicity Z = 10“"^ (that is [Fe/H] ^ 


-2.3) with their population-synthesis models, adopting the solar- 
neighbourhood IMF proposed by |Kroupa et al.| ( 19931, and vary¬ 
ing a set of uncertain physical parameters related to the mecha¬ 
nism of mass transfer, nucleosynthesis, and mixing processes. 
They typically find CEMP fractions between 2% and 4%, while 
the observed values found by different authors at [Fe/H] < -2 
vary between 9% and 25% ( Marsteller et al.|2005[[Frebel et al.[ 


[2066[ [Lucatello et al.[[2006 [Lee et al.[[2013| l. However, their 
AGB evolution models are based on detailed models that do not 
allow for third dredge-up in AGB stars of mass below approx¬ 
imately 1.2 Mq. To increase the CEMP/VMP fraction, [Izzard 
et al. ( 2009| l and Abate et al. ( 2013| l artificially allow for third 


dredge-up in AGB stars down to mass 0.8 Mq in their models (as 
proposed by Izzard et al. 2004 using independent arguments), 
and show that in this way it is possible to raise the modelled 
CEMP/VMP ratio to > 10%. 

The most recent detailed models of AGB evolution and nu¬ 
cleosynthesis of Stancliffe & Glebbeek| POO’S] !, Karakas ( 2010| l 
and[Lugaro et al.[([2012|l find third dredge-up in stars of mass 


down to 0.9 Mq at metallicity Z = 10 In the models of 
Karakas (20101 and Lugaro et al. (2012|l, the nucleosynthesis 


products of stars are computed on a grid of 16 initial masses in 
the range [0.9, 6] Mq, taking in to account 320 isotop es from 'H 
up to ^'°Po. In our recent work (Abate et al. 2015a|b Paper I and 
II hereafter) we include the results of these detailed calculations 
in our model of binary evolution and nucleosynthesis, and study 
in detail the chemical compositions for a sample of 67 CEMP-s 
stars. We compare the surface abundances of each CEMP-s star 
with a grid of model binary stars, and determine the best fit to the 
observed abundances. From this comparison we conclude that 
our models reproduce reasonably well the chemical properties of 
CEMP-s stars with [Eu/Fe] < +1 (for brevity, “r-normal CEMP- 
s” stars hereafter), but not CEMP-s/r stars, possibly indicating 
different nucleosynthesis histories for the two categories. Also, 
our models typically predict orbital periods longer than those 
observed, suggesting that wind mass transfer should be more ef¬ 
ficient at close separations. 

In this paper we extend our recent work on CEMP stars, 
and analyse the properties of the population of CEMP stars at 
[Fe/H] ai -2.3 as a whole, through a comparison with synthetic 
populations. Our purpose is to answer the following questions: 

1. Is it possible to reproduce the observed CEMP/VMP ra¬ 
tio with our model of binary population synthesis that in¬ 
cludes the latest detailed AGB models of [Karak^ ( |2010| l and 
Lugaro et al. (20121? 


^ The original nomenclature for these stars is “CEMP-r/j”; we em¬ 
ploy CEMP-5/r for consistency with Papers I and II. 


2. Can we constrain our model of the wind mass-transfer pro¬ 
cess, and in particular the wind-accretion efficiency and 
angular-momentum loss? 

3. Under which conditions does our model reproduce the ob¬ 
served abundance distributions of elements associated with 
AGB nucleosynthesis? 

The paper is organised as follows. In Sect. [^ we describe 
the most important characteristics of our model and the selec¬ 
tion criteria that we adopt to compare the results of our synthetic 
population with the observations. In Sect. [^ we present the re¬ 
sults of our models, which are further discussed in Sect. [4[ Our 
conclusions are presented in Sect. [^ 

2. Models 

In this work we use the binary population synthesis code 
binary_c/nucsyn described by |Izzard et al.[ ( |2004|[2006[[2009| l 
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to simulate populations of stars at low metallicity, and compare 
them with a sample of observed stars. In this section we sum¬ 
marise the physical parameters of our model (Sect. |2.1|i, th e as¬ 
sumptions made in our grid of binary models (Sect. |2.2) i, and 
describe the method used to select the simulated stars according 
to their luminosity (Sect.[23]l. 


2.1. Input physics 

In our default model set A we adopt the same input physics as in 
model set A of Papers I and II. We briefly list our most important 
assumptions below. Table[^provides an overview of our alterna¬ 
tive model sets, in which we varied some of these assumptions. 


The Reimers ( 1975| l equation multiplied by a factor of r] - 
0.5 is used to compute the wind mass-loss rate up to the 
AGB. The prescription of Vassiliadis & Wood] ( 1993) l de¬ 
scribes mass loss during the AGB phase, with minimum 
and maximum wind velocities = 5kms“' and Vw = 
15 kms^', respectively. 

We adopt the approximation of a spherically-symmetric 
wind ( Abate et al.|2013 Eq. 4) to compute the angular mo¬ 
mentum carried away by the ejected material, and calcu¬ 
late the wind-accretion rate according to a wind Roche-lobe 


overflow (WRLOF) model described by Eq. (9) of Abate 
et al.| (|2OT3 i. Our model sets B and C are as in Paper II. 


In both model sets the expelled material carries away a mul 
tiple 7 = 2 of the average specific orbital angular momen¬ 
tum of the binary system. In model set B we adopt an en¬ 
hanced Bondi-Hoyle-Lyttleton (BHL) model of wind accre¬ 
tion computed with Eq. ( 6 ) of |Boffin & Joriss^ ( 1988] l and 
Q'bhl = 10, whereas in model set C our default WRLOE pre¬ 
scription is adopted (Paper II). 

The algorithms that we use to compute the nucleosynthe¬ 
sis of the star through the first and second dredge-ups are 
based on the results of [Karakas et al.| (|2002[) and [Karakas 


& Lattanzio 


et al. (2004 


(|2007 1; we refer the interested reader to |Izzard 


2006 2009 [) for a thorough description. 


The amount of material dredged up during the AGB phase 
is compute d with algorithm s tuned to reproduce the detailed 
models of Karakas (2010|l at metallicity Z = 10 as de¬ 
scribed in Paper I. The chemical composition of the inter¬ 
shell region of the AGB star is tabulated as a function of 
the mass of the star at the beginning of the AGB phase, the 
evolution along the AGB, and the mass of the partial mix¬ 
ing zone, MpMZi ^ free parameter of our model. In the de¬ 
tailed AGB models of [Karakas ( 2010| l, the partial mixing 
zone (PMZ) is the site in which protons from the envelope 
are partially mixed with material of the intershell region. The 
protons react with the nuclei and form a pocket, 
in which free neutrons are produced and become available 
for the i-process. We refer the reader to Karakas 1^010) 1 and 
Lugaro et al.|(2012]l for a detailed explanation of the numer¬ 


ical treatment of the PMZ, and to Paper I for a description of 
our implementation in our population-synthesis code. 

We assume efficient thermohaline mixing, that is, the ac¬ 
creted material is instantaneously mixed with the accret¬ 
ing star. The calculations of Stancliffe et al. ( 2007| l suggest 
this approximation is reasonable in many cases, even though 
[Stancliffe & Glebbeek] ( |2008| l and jStancliff^ ( |2009| ) show 
that other processes, such as gravitational settling, may in 
some circumstances reduce the effect of thermohaline mix¬ 
ing. To account for the possibility of inefficient thermohaline 


Table 1. Physical parameters adopted in our binary population 
models. 


model set 

Physical parameters 


(differences from model set A) 


A 

Al 

A4 

B 


D 

Ikoi 

Il05 

Qi 

Qp4 

Sdm 

S3 

T8 

T12 


Default 

[C/Fe] > 1.0 

As Al, logjQ(g/cms“^) < 4.0, no luminosity selection 

Angular-momentum loss: AJ/J = 2(AM/M) 
Wind-accretion efficiency: BHL, asHL = 10 
(as model B of Paper I and II) 

Angular-momentum loss: AJ/J = 2 (AM/M) 
Wind-accretion efficiency: WRLOF 
(as model C of Paper II) 

No thermohaline mi xing 


Multiple-part power-law IMF (Kroupa 2001 1 , 
forM > IMq , IMF slo pe is less steep ttian d etault. 
Log-normal IMF i Lucatello et al. 2005 i, 

PlogM ~ 0.79, (TiQgj^ — 1.18. 

^(q) oc q 
(/>(q) oc q-°-* 

_ Log-normal initial-p eriod distribution 

Duquennoy & Mayor] 1991 1 , yUiog/> = 4.8, criogp - 2.3. 
As Sdm with piogp = 3.0 
tmin “ 8 Gyr 
tmin = 12 Gyr 


mixing, in model set D the accreted material remains on the 
stellar surface until mixed in by convection. 

In our models we assum e the metallic i ty of t he detailed nu¬ 
cleosynthesis models by Lugaro et al. ( 2012| l, i.e. Z = 10 
We adopt the same initial composition as Papers I and II 


based on the chemical evolution models of Kobayashi et al. 
(2011 1 for the isotopes up to ^®Ge. For heavier isotopes we 


assum e the solar abundance distribution of Asplund et al. 
(|2009ll scaled down to Z = 10 


2.2. Population synthesis 

Each of our simulated populations consists of N binary stars uni¬ 
formly distributed in the logjo Mi -M 2 -logjQ a parameter space, 
where Mi 2 are the masses of the primary and the secondary star, 
respectively, and a is the orbital separation. Seven different val¬ 
ues of the PMZ mass are taken into account. The adopted grid 
resolution is A = Ami x A^m2 x A^ x Apmz, where Ami = 100, 
Am2 = 32, Na - 80 and Apmz = 7. A finer grid with higher 
resolution does not change our results. 

The initial parameters vary in the following ranges: 


Mi l is in the interval [0.5, 8.0] Mq. Stars of mass below 
0.5 Mq are not expected to be visible with the current mag¬ 
nitude limitation of the observational surveys, as explained 
in Sect. 2.3 Massive stars that end their evolution as super¬ 
novae are not considered. 

M 2 ,i is uniformly spaced in [0.1, 0.9] Mq. More massive stars 
have already evolved to a white dwarf after 10 Gyr, thus they 
are excluded because they do not contribute to the fraction 
of CEMP stars. By definition, initially M 24 < Mi i, therefore 
the initial mass ratio qi - M 2 ,i/Mi,i is always 0 < q'l < 1 . 

We consider circular orbits with ai between 50 and 5x IO^Rq. 
Binary stars outside the considered range do not become 
CEMP stars in our models, either because they are too close 
and they merge when the primary star becomes a red gi- 
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ant, or because they are too wide to interact. However, bi¬ 
nary stars at closer separations are observed. For exam¬ 
ple, |Kouwenhoven et al ( 2007| l analyse the binary popu¬ 
lation of the stellar association Scorpius OB2 in the range 
[5, 5 X 10^] Rq. We assume that all stars are formed in binary 
systems with separations in this range. When calculating the 
CEMP fraction, we take into account that our grid investi¬ 
gates a fraction of the total population of binary stars. For 
example, with a flat distribution in logjg fli our grid covers 
83% of the observed range, thus the model CEMP fraction 
is multiplied by this factor. Eor different ai distributions we 
compute the corresponding corrections. Binary systems with 
wider separations than 5 x 10^ Rq and single stars could in 
principle be taken into account by decreasing the binary frac¬ 
tion. 

When the stellar mass is M* < 3 Mq we adopt the following 
values of Mpmz: 0, 2 x 10-\ 5 x IQ-^^, 6.6 x 10-^ 1 x 10-^ 
2 X 10“^, 4 X 10“^ Mq. Mpmz is always zero for M, >3 Mq. 


Pollowing |Izzard et al. ( 2009| l and Abate et al. (2013 1 , stars 
are counted according to the sum 


A^l,max A^2,max <3n,ax ^max 

Wtype ^ X Xi Xi 6M2 5a 6t , (1) 

^2,min ^min ^min 


where 

- S is the star-formation rate, assumed to be constant for ages 
between fmin and fmax- 

- fmin = 10 Gyr corresponds to the approximate minimum age 
of the Galactic halo, while f^ax = 13.7 Gyr is the age of the 
Universe; 6t is the timestep. 

- The volume unit in the parameter space is 5Mi ■ 6M 2 ■ 5a. 

- <Jtype = 1 when the star belongs to a specific type, and 
zero otherwise. Model stars are classified as follows; VMP 
stars are all stars older than tmin that have not become white 
dwarfs; CEMP stars are VMP stars with [C/Ee] > h- 0.7 ; 
CEMP-i stars are CEMP stars with [Ba/Ee] > -1-0.5. 

- 'Em!,M 2 ,a is the inital distribution of Mi, M 2 , and a. We as¬ 
sume to be separable. 


T' = T'(Mi,M2,fl) = (^(Mi) 4>iM2)x(a) ■ 


( 2 ) 


In our default model set A the primary mass distribution 
tl/(Mi i) is the IME proposed by |Kroupa et al.| (|1993|l, the 
secondary mass distribution 0 (M 2 ,i) is flat m qi, and the sep¬ 
aration distribution is flat in log a. 


Some of the assumed parameters in our model are uncertain 
and not well-constrained by observational data, in particular: the 
IME of the primary star, the binary fraction, the distribution of 
the mass ratios and separations in binary systems, and the age 
of the halo. To determine how each of these uncertainties affects 
our results we modify our default model set A by varying one 
parameter at a time. Table[T]lists some of the model sets that we 
tested, and their differences with respect to our default model 
set A. 


2.3. Selection criteria 

In the papers by |Izzard et al.| ( 2009[ ) and [Abate et al.| ( |2013| l, 
the synthetic stars are selected from the simulations when their 
surface gravities are below the threshold logio(g/cms“^) = 4.0. 
This criterion essentially restricts the analysis to giant stars, and 


is based on the implicit assumption that all giant stars of the halo 
are visible. Here we replace this criterion with a more realistic 
selection, based on the apparent magnitudes of the stars, includ¬ 
ing in the analysis the main-sequence stars that are sufficiently 
luminous to be visible. 

To calculate the number of stars that are visible as a func¬ 
tion of their luminosity, N{L), we follow the method adopted by 
van Orrschot et al. ( 2014| l. The main steps are summarised as 
follows. 


1. Our model computes the luminosity, L, the effective tem¬ 
perature, Teff, and surface gravity, logjg g, of each simulated 
star. The absolute visual magnitude. My, is calculated from 


the luminosity as described by Torres (2010 Eq. 10). The 
bolometric correction is computed as a function of T^ff and 
login g of the star, adopting the values published by [Bessell 
jetaLj ( [19981 Table 1). 

2. The apparent V magnitude of a star is computed as a function 
of the distance from the Sun, d, and the Galactic latitude, 
b. Similarly to van Orrschot et al. ( |2014| l, the extinction of 
stellar light is taken into account according to the Galactic- 
dust distribution prop osed by [Toonen & Nelemansj ( |2013| l 
and the prescription of Sandage 1 19721. For ^4 0, one has 


V = My + 5 (logio(^^) - 1) + Ay(oo) tanh 


jdsinb 

\ Zh 


1 (3) 


where Zk - 120 pc is the scale height of Galactic dust, and 
Ay(oo) is: 


Ay(oo) = 


0.165(tan50°-tanZ7)cscfe if |fe| < 50° 


0 


if\b\ > 50°. 


3. For each simulated star of luminosity L and corresponding 
absolute magnitude My, we compute the volume in which 
the star has apparent magnitude < V < Vmax- In our 
simulations, Vmin and Umax are set consistently with the de¬ 
tection limits of the observational data. We compare our syn¬ 
thetic populations with two different data sets, which have 
different magnitude ranges: 

i. The observations of the SDSS/SEGUE survey (for 
brevity, SEGUE), which include the effective tempera¬ 
tures, surface gravities, and carbon abundances for about 
13,000 very metal-poor stars with-2.5 < [Ee/H] < -2.0 
observed at medium resolution ( |Lee et al.|2013] l. The for¬ 
mal g-magnitude range of the SEGUE survey is 14 < 
g < 20.3; however, Yanny et al. (2009 1 state that reliable 
atmospheric parameter measurements are possible only 
within the range 14 < g < 19. We convert the limits in 
the g-band to 13.5 < V < 18.5 with the empirical photo¬ 
metric relation ([Windhorst et al.|1991|). 


y = g - 0.03 - 0.42 (g - r). 


(4) 


where we adopt (g - r)niax = 1 -2, as derived from Eig. 7 
of [Lee et al.[ ( |joT^ . 

A sample of 378 very metal-poor stars with -2.8 < 
[Ee/H] < — 1 . 8 , both carbon-normal and carbon- 

enhanced, mostly based on the SAGA database (|Suda[ 


et al. 2008[ 


2011 last updates in January, 2015). This 


is the sample adopted in Paper II, and includes high- 
resolution measurements of the surface abundances of 
many elements. The apparent V magnitudes of the stars 
in this sample vary between Vinin = 6 and Fmax = 16.5, 
hence we adopt this range in our simulations. 
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4. To determine N{L, logjQ g), that is, the number of visible 
stars with luminosity L, temperature Teff, and surface grav¬ 
ity logjo g, we integrate the stellar density distribution of the 
Galactic halo over the volume determined in the previous 
step. We assume a maximum distance of d - 10^ pc. The 
stellar density distribution in the halo is parameterised with 
the equation of an oblate spheroid in the reference frame of 
the Galactic centre as prescribed by Helmi|([20081l, 


Table 2. CEMP/VMP ratio, Fq, calculated with model sets Al 


and A4 and with set Cq of 

Abate et al. 

(2013 

1 . The errors con- 

vey only Poisson statistics 

The SEGU 

b value ot I^c the 


uncertainty are taken from|Lee et al.|([2()13[ Tab. 4) for stars in 


the range -2.5 < [Fe/H] < -2. Stars from the SAGA database 
are selected with iron abundance in the range [- 2 . 8 , - 1 . 8 ] to in¬ 
crease the number statistics. In all cases CEMP stars are dehned 
as stars with carbon abundance [C/Fe] > -i-l. 


p(x, y, z) oc Tq" -b -I- (azf )" , (5) 


where Tq = 8 kpc is the distance of the Sun from the Galactic 
centre (e.g. |Moni Bidin et al.|2012] l, a is the minor-to-major 
axis ratio, and n is the exponent of the density prohle. For a 


and n we adopt the best-fit values determined by Juric et al. 
( |2008| l in their three-dimensional density map of the Galactic 
halo: a - 0.64 and n — -2.8. 


5. To speed up the calculations we initially calculate 
A(L, Teff, logjog) for a grid of luminosities, temperatures, 
and gravities. The intervals of these parameters in the 
grid are chosen to reproduce the ranges of variation of 
these parameters in the binary stars of our simulations 
(namely: [-3, 5], [3000, 8000] and [-1, 6 ] for log[o(L/Lo), 
Teff/ K and logio(g/cms“^), respectively). The values of 
N{L, Teff, log [0 g) are stored in two tables, one for each mag¬ 
nitude range considered. 

6 . Each star in our simulation is counted by multiplying the 
value Utype computed with Eq. (fH by N(L, T^s, logjg g). The 
value of A(L, 7’eff,logjQg) is determined by interpolating 
L, Teff, and logjQ g within our tables. This step is repeated for 
all the synthetic populations generated with different model 
sets. 


model sets and data 

Fc (%) 


Al 

8.05 ± 0.02 


A4 

7.91 ±0.02 

Cqi 

Abate et al. 2013^ 

2.38 ± 0.04 


SEGUE 

6.1 ± 1.0 

SAGA database 

23.5 ± 2.9 



Finally, we note that in our models most horizontal-branch 
stars have large effective temperatures, up to 10,000 K. In the 
range -2.8 < [Fe/H] < -1.8, the SEGUE and SAGA databases 
contain only a few horizontal-branch stars hotter than 6,000 K, 
possibly because hotter stars are selected against. In our simu¬ 
lations these stars are typically not very numerous because the 
time spent on the horizontal branch is short compared to the 
total stellar lifetime. However, including them in our analysis 
modifies the model distribution of surface gravities between ap¬ 
proximately 2.5 dex and 3.5 dex, as noted also by [Izzard et al. 
( |2009| l. For this reason we exclude from our simulations all stars 
with Teff > 6,000 K that have evolved further than the turnoff. 
This selection has a negligible effect on the computed fraction 
of CEMP stars. 


3. Results 


3.1. Comparison with the resuits of \Abate et at.\(20i:^ . 


We now compare the CEMP-to-VMP fraction, Fq, computed 


with our default model set A and model set Cq of [Abate 


et al. ( 2013| l. These two model sets adopt the same assump¬ 
tions about the distributions of initial parameters, the wind- 
accretion efficiency, and the mechanism of angular-momentum 
loss. However, the selection criteria used in these model sets 
are different (Sect. |2.3| l. In model set A4 the same dehni- 
tions are used as in model set Cq to select VMP stars, namely 
logio(g/cm s“^) < 4.0 and f* > 10 Gyr, and the criterion 
[C/Fe] > -bl is used to select CEMP stars. Because model set 
Cq was computed assuming a range of initial separations differ¬ 
ent from the present work, namely between 3 Rq and 10^ Rq, the 


Fig. 1. Distribution of Mi i computed with model set A4 (solid, 
blue line) and model set Cq of Abate et al. (2013 1 (dashed, red 
line). The bins are equally spaced and the width in Mi i/Mq of 
each bin is 0.1. The y-axis indicates the number of CEMP stars 
in each bin, arbitrarily rescaled. 


CEMP/VMP ratio of model Cq has been recalculated to take into 
account this difference. 

With these assumptions, the only difference between model 
sets A4 and Cq is in the adopted set of AGB evolution models. 
The AGB model of Abate et al. (20131 was tuned to reproduce 
the evolution of AGB stars as predicted by the detailed mod¬ 
els of jKarakas et al.| ( [2002| l and jKarakas & Lattanzio|p007| l. In 
these models, AGB stars of mass below approximately 1.2 Mq 
did not experience third dredge-up, and consequently did not 
contribute to the formation of CEMP stars. On the other hand, 
jStancliffe & Glebbeekj (|2008|l and the recent work of [Karakasj 
020101 1 , on which our AGB model is based, show that AGB stars 
of mass down to 0.9 Mq undergo third dredge-up. Consequently, 
in our model set A4 CEMP stars are formed in a wider range 
of Miq. The distributions of Mi i computed with model sets A4 
and Cq are compared in Fig. Model set A4 peaks for Mi i in 
the range [0.9, 1.1] Mq, in which set Cq does not produce CEMP 
stars. With the adopted solar-neighbourhood IMF low-mass stars 
are much more commo, consequently model set A4 predicts a 
CEMP fraction of 7.9%, i.e. a factor 3.3 larger than set Cq, as 
shown in Table |2] 

Model set Al differs from model set A4 in that the model 
stars are selected according to their luminosity, as described in 
Sect. |2.3[ With model set Al we hnd a CEMP/VMP ratio of 
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8.0%. Comparison with the result of model set A4 indicates that 
the selection criterion has a small impact on the total fraction 
of CEMP stars because in our synthetic population the fraction 
of CEMP stars is very similar among stars above and below 
the threshold of logiQ(g/cms“^) = 4. If we reduce the mini¬ 
mum abundance of carbon in the definition of CEMP stars, and 
adopt [C/Ee] > h- 0.7 (our default assumption in this paper), 
the CEMP/VMP ratio increases to 8.5% (model set A, Table 
[^. This occurs because our model produces very few stars with 
-1-0.7 < [C/Ee] < -hI.O, whereas 94% of the model CEMP stars 
have [C/Ee] > -i-l, and the average carbon abundance in the syn¬ 
thetic CEMP population is [C/Ee] = -1-1.9. 

Table compares the model predictions with the observed 
CEMP/VMP ratio derived from the SEGUE and SAGA data 
sets. The CEMP fraction among SEGUE stars is taken from 
Table 4 of [Lee et ah] ( |2013| l for stars with iron abundance in 
the range -2.5 < [Ee/H] < -2 and with [C/Ee] > -i-l. A more 
detailed comparison between the results of our model sets and 
the SEGUE sample is performed in the next section. The CEMP 
fraction in the SAGA sample is approximately 24%, almost four 
times larger than the value determined for SEGUE stars. We dis¬ 
cuss this discrepancy in Sect. 


3.2. The fraction of CEMP stars 


We now compare the CEMP fraction observed in the SEGUE 
sample with the results of our simulations. We select stars ac¬ 
cording to the method explained in Sect. |2.3[ with Umin = 13.5, 
Umax = 18.5 and a carbon abundance in CEMP stars [C/Ee] > 
-1-0.7. 

Table|3summarises the results obtained with different model 
sets. The observed fractions and relative uncertainties are com¬ 
puted from the data reported by [Lee et ah] (|2013| Table 6) for 


stars in the range -2.5 < 


assumptions of Lee et al. 


Ee/H] < -2. Eor consistency with the 


(20131, a minimum observational un¬ 


certainty of 1% is adopted, and stars are defined as giant, turnoff, 
and dwarf if their surface gravities are log[Q(g/cms“^) < 3.5, 
3.5 < log[o(g/cms“^) < 4.2 and log[o(g/cms“^) > 4.2, respec¬ 
tively. The CEMP/VMP ratios computed among giants, turnoff, 
and dwarf stars are indicated as Ec,gia, Fcjo, Ec,dwa (columns 
3 - 5). The proportions of giant, turnoff and dwarf stars among 
CEMP stars are indicated as Ugia.cEMP, Eto.cemp, Udwa.CEMP, re¬ 
spectively (columns 6 - 8). The quoted uncertainties in the mod¬ 
els convey only Poisson statistics. In this section we compare the 
observed and model fractions calculated with our default model 
set A, while the comparison of the results with different model 
sets is the subject of the next section. 

Our total CEMP/VMP fractions underestimate those ob¬ 
served by a factor of 1.2 - 1.4. We note a large difference be¬ 
tween the observed and modelled proportions of the three evo¬ 
lutionary stages among CEMP stars. In our models almost half 
of the CEMP stars are on the main sequence. By contrast, in the 
observed sample the CEMP dwarf stars count for approximately 
1 .6%, while giants and turnoff stars count for approximately 
50% each. This discrepancy is also evident from the compari¬ 
son between the logjy g distribution of the observed CEMP stars 
of the SEGUE sample and the results of model set A (Eig. |^. 
The positions of the peaks of the two distributions differ by ap¬ 
proximately 0.5 dex. If we artihcially decrease the logjQ g of our 
model CEMP stars by 0.5 dex to mimic a hypothetical system¬ 
atic offset we obtain fractions of giant, turnoff, and dwarf stars 
among CEMP stars of 35.4%, 60.5% and 4.1%, respectively. 
However, with such a systematic offset we hnd that approxi¬ 
mately 3% of the modelled CEMP stars have log[o(g/cms“^) < 



Fig. 2. Distribution of log[g(g/cms as derived from the 
SEGUE sample of CEMP stars with -2.5 < [Ee/H] < -2 
(histogram with Poisson errors) and computed with our default 
model set A (solid line). The 0.2 dex bins are equally spaced in 
log[g(g/cm s“^). In this and the similar plots that follow the area 
under the graph (which is the total number of stars) is normalised 
such that it is the same for both the observations and our model 
stars. 


1, while far fewer SEGUE stars are observed at such low gravity 
QLee et al.|20131 l. 

An offset in the model-vs.-observed surface gravities can 
also be identified in Eig.|^ in which the temperature-gravity dis¬ 
tribution of our synthetic population (red-scaled distribution) is 
compared to SEGUE stars with iron abundances in the range 
-2.4 < [Ee/H] < —2.2. The surface gravity of the main se¬ 
quence and of the turnoff point in the synthetic population ap¬ 
pear to be shifted by several times 0.1 dex with respect to the 
observed sample. The reason for such a discrepancy is currently 
unclear. One possible explanation is that in our simulations we 
overestimate the upper limit in V magnitude and, as a conse¬ 
quence, in our selection procedure we take into account stars 
that in reality are too faint to be detected. However, Table 
shows that, even if Umax is reduced by 2 mag, the proportion 
of dwarfs among CEMP stars decreases only by a factor of two. 
To reproduce the observed dwarf fraction a far lower Umax is 
necessary, which is inconsistent with the characteristics of the 
SEGUE survey. An alternative explanation of this discrepancy 
is that the surface gravities determined from the SEGUE spectra 
are systematically underestimated. However, the algorithms to 
determine the surface parameters of SEGUE stars are calibrated 
on a set of stars with available high-resolution spectroscopy, and 
systematic errors that shift logipg by mor e than -0.3 dex are 
considered highly unlikely ( Lee et al.|2013 1. 

Eigurej^ shows that there is also a discrepancy between the 
effective temperatures of the modelled and observed stars. To 
reliably estimate [C/Ee], [Lee et al.| ( |2013 l restrict the SEGUE 
sample to stars with 4,400 K < T^g < 6,700 K because in this 
range their technique is robust. However, in our models most of 
the turnoff stars have Teff > 6,700 K, up to « 8,000 K, regardless 
of the predicted abundance of carbon, whereas in the SEGUE 
sample there are very few turnoff stars with T^g > 6,700 K 
(which are not plotted in Eig. [^. Such a discrepancy is at least 
partly explained by an inaccuracy in our model temperatures. 
Our code calculates the effective temperature from the luminos¬ 
ity and the radius of the star, which are computed with htting for- 
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Fig. 3. Feff - logiQ^ probability distribution of all stars in our 
simulation computed with model set A. Darker colours in¬ 
dicate regions of higher probability. The grey solid and ma¬ 
genta dotted lines show the isochrones computed with detailed 
stellar-evolution models at 10 Gyr and 13.7 Gyr, respectively. 
The isochrones computed for the same ages with a detailed 
model that includes gravitational settling are shown as light-blue 
dashed and blue dot-dashed lines, respectively. The black dotted 
lines indicate the ranges of and log mg in which |Lee et al. 
( |2013| l select the SEGUE stars. Black circles indicate all SEGUE 
stars with -2.4 < [Ee/H] < -2.2. 


mulae based on the detailed stellar models of |Pols et al.| ( |1998l l. 
Eor a 0.8 Mq star at the turnoff point these fitting formulae re¬ 
produce the luminosity and radius of the detailed model within 
a factor of 1.15 and 0.98, respectively. However, because with 
binary_c the luminosity is higher and the radius is smaller than 
in the detailed model, the combination of these two factors gives 
a higher effective temperature than in the detailed model. Eig. 
shows two isochrones, for 10 Gyr and 13.7 Gyr (grey solid and 
magenta dotted lines), which were constructed using the same 
version of the STARS code (|Eggleton||197H subsequently up¬ 
dated by many aut hors, e.g. |Pols et al.|1995| as adopted in the 
detailed models of |Pols et al.| ( |1998| l. These isochrones bracket 
the age range of our synthetic population. The offset between 
our population-synthesis results and these isochrones yields an 
estimate of the error in temperature because of inaccuracies in 
the fitting formulae. 

The discrepancy between modelled and observed tempera¬ 
tures decreases if efficient gravitational settling is included in 
the models. Two isochrones including settling are shown in 
Pig. 1^ (light-blue dashed and blue dot-dashed lines, for 10 Gyr 
and 13.7 Gyr, respectively), which were constructed using a re¬ 
cent version of the STARS code ( |Stancliffe & Glebbeek||20()8] 
IStancliffe & Eldridge|2009] l. Settling decreases the turnoff tem¬ 
perature of the 10 Gyr isochrone by approximately 300 K, that 
is, the same shift caused by a 3 Gyr difference in age. The 
13.7 Gyr isochrone with settling has a turnoff temperature of 


about 6400 K, which is not much hotter than the average turnoff 
temperature of the observed stars (=» 6000 K). Although the ef¬ 
fects of gravitational settling on the structure and composition of 
metal-poor stars are not fully explored ( [Richard et al.|2002a ^ 
Stancliffe & Glebbeek|2008 1, these results suggest that the dis¬ 


crepancy in effective temperature can be reduced by including 
settling in our models. 

The effective temperature of giant stars is underestimated 
both by the population synthesis models and by the isochrones. 
This discrepancy can also be resolved by calibrating the mixing- 
length parameter in the detailed models (as done for example by 
Spada et al.|2013 in the Yonsei-Yale isochrone sets), and it does 


not affect any of the results presented in this work. 

In the observed sample the CEMP/VMP ratio is found to in¬ 
crease signihcantly with luminosity, that is, Fq is smaller among 
dwarf stars and larger among giants, although this relation is re¬ 


versed at lower metallicity (Pe/H < -3, Lee et al. 


2013 I. In 


our models we find smaller variations of the CEMP/VMP ratio 
among the different classes of stars. The smallest CEMP/VMP 
ratio is calculated among dwarf stars and the largest among 
turnoff stars. Turnoff and giant CEMP stars in the models have 
very similar mass distributions: the CEMP fraction is smaller 
among the giants as an effect of the hrst dredge-up, which 
depletes part of the carbon at the surface. Consequently, a 
turnoff CEMP star with carbon abundance close to the threshold 
[C/Pe] = h- 0.7 has [C/Pe] < h- 0.7 after the dredge-up. This ef¬ 
fect is small because in our model thermohaline mixing already 
dilutes the accreted material throughout the entire star. If we in¬ 
hibit thermohaline mixing in our models the difference between 
^c.gia Fq jq increases (Table model set D). 

The reason for the discrepancy between the predictions of 
our models and the observations is still unclear. One possible 
explanation is related to the method adopted by Lee et al. ( |20131 l 
to compute Fq. Stars for which only an upper limit to [C/Pe] 
is available or stars with a poor carbon measurement (i.e. if the 
correlation coefhcient between the observed and synthetic spec¬ 
trum is lower than 0.7) are counted among carbon-normal very 
metal-poor stars. Carbon lines are more difficult to detect in the 
spectra of warm stars, which are therefore clearly identified as 
CEMP stars only if carbon is very enhanced. As a consequence, 
the CEMP/VMP ratio computed for turnoff stars is somewhat bi¬ 
ased towards lower values, although this effect is probably small 
( Lee et al.|2014| ). This bias is not expected to be signihcant for 
the observed dwarf stars, which are typically cooler than turnoff 
stars. The fact that the observed Ec.dwa is approximately a factor 
of ten lower than Fc.to and Fc.gia suggests some effect prevents 
the detection of dwarf CEMP stars. It is still unclear whether 
this is a physical effect, which should be stronger for carbon- 
rich dwarfs than for carbon-normal ones, or a bias introduced 
during the spectral analysis. 


3.3. CEMP fractions computed adopting different sets of 
initiai parameters 

The characteristics of our synthetic CEMP populations depend 
on our assumptions about the mass-transfer process, the extent to 
which the accreted material is diluted, and the properties of the 
stellar population, such as its age, and the initial distributions of 
masses and separations. 

The differences in the CEMP/VMP ratios between model 
sets A, B, and C are mainly due to the range of periods at which 
the secondary star accretes material from the AGB primary star, 
as shown in Pig. Model sets B and C have a larger fraction 
of CEMP stars than set A because of the more efficient mecha- 
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Table 3. CEMP star fraction computed from the SEGUE data in the range -2.5 < [Ee/H] < 2, and computed with our different 
model sets. The range of V magnitudes adopted in the models is [13.5, 18.5]. Fq, Fc,gm, Fqjo, Ec.dwa indicate, respectively, the 
CEMP/VMP ratio computed among all/giant/turnoff/dwarf stars. The fraction of giants, turnoff stars and dwarf stars computed 
among CEMP stars is indicated respectively as Ugia.cEMP, Eto.cemp, Udwa.CEMP (the latter fractions add up to 100%). The errors in 
the models convey only Poisson statistics. A minimum observational uncertainty of 1% is adopted (cf. Lee et al.|2013| l. CEMP stars 
have carbon abundance [C/Ee] > +0.7. 


model sets and data 

Fc (%) 

Cc,s.a(%) 

F C.TO (%) 

Cc.dwa (%) 

F gia.CEMP (%) 

F TO.CEMP (%) 

F dwa.CEMP (%) 

SEGUE 

11.46+1.00 

15.79+ 1.00 

10.38+ 1.00 

1.67+1.00 

51.42 + 2.30 

47.02 + 2.16 

1.56+1.00 

A 

8.50 ± 0.02 

9.59 + 0.03 

10.73 + 0.04 

7.32 + 0.02 

31.95 + 0.10 

19.80 + 0.08 

48.25 + 0.13 

B 

8.96 ± 0.02 

10.30 + 0.03 

11.01+0.04 

7.72 + 0.02 

32.47 + 0.10 

19.23 + 0.07 

48.30 + 0.13 

C 

9.48 ± 0.02 

10.67 + 0.03 

11.78 + 0.04 

8.24 + 0.02 

31.69 + 0.10 

19.56 + 0.07 

48.75 + 0.13 

D 

12.17 + 0.02 

9.94 + 0.03 

15.21+0.05 

12.45 + 0.02 

23.12 + 0.07 

19.60 + 0.06 

57.28 + 0.12 

Ikoi 

8.81+0.01 

9.98 + 0.03 

11.12 + 0.04 

7.58 + 0.02 

31.99 + 0.10 

19.75 + 0.08 

48.27 + 0.13 

Il05 

16.47 + 0.02 

17.23 + 0.04 

18.31+0.05 

15.40 + 0.03 

32.23 + 0.08 

19.27 + 0.06 

48.50 + 0.10 

Qp4 

7.87 + 0.01 

8.69 + 0.03 

9.85 + 0.04 

6.90 + 0.02 

31.37 + 0.11 

19.45 + 0.08 

49.18 + 0.14 

Qi 

11.37 + 0.02 

12.18 + 0.03 

13.95 + 0.05 

10.11+0.02 

32.27 + 0.10 

20.28 + 0.07 

47.45 + 0.13 

Sdm 

11.16 + 0.02 

12.45 + 0.03 

14.06 + 0.04 

9.67 + 0.02 

31.98 + 0.08 

19.72 + 0.06 

48.29 + 0.11 

S3 

9.41+0.02 

10.72 + 0.03 

11.82 + 0.04 

8.08 + 0.02 

31.71+0.09 

20.05 + 0.07 

48.24 + 0.12 

T8 

9.89 + 0.01 

12.14 + 0.03 

10.11+0.03 

8.78 + 0.02 

30.91+0.09 

20.22 + 0.07 

48.86 + 0.12 

T12 

8.14 + 0.02 

9.72 + 0.04 

13.60 + 0.08 

6.43 + 0.02 

31.52 + 0.16 

19.70 + 0.12 

48.78 + 0.21 


nism of angular-momentum loss which allows binary systems to 
interact at longer initial periods. The largest range of initial peri¬ 
ods is accessible to model set C because it combines the model 
of efficient angular-momentum loss with the WRLOE model of 
wind accretion. Similar results are found by Abate et al. ( 2013| l. 
In binary systems that are initially very wide (P; > 10° days) 
only relatively massive primary stars (Mi^ > 1.5 Mq) transfer 
significant amounts of mass. Such stars are rare according to the 
solar-neighbourhood IME adopted in our models. Consequently 
the increase in Fq in models B and C is rather modest. 

The relative proportion of giant, turnoff, and dwarf stars de¬ 
pends essentially on the mass distribution of the CEMP stars. All 
model sets produce the majority of CEMP stars with masses af¬ 
ter accretion between 0.75 Mq and 0.9 Mq. However, the initial 
mass distribution of the CEMP stars and the distribution of the 
accreted masses differ between model sets. In wide binary sys¬ 
tems {P > 10^ days) our WRLOE model typically predicts larger 
accretion efficiencies than our enhanced Bondi-Hoyle-Lyttleton 
model. Hence, at large separations secondary stars on average 
accrete more mass with model sets A and C than does model set 
B, and are initially less massive. As a consequence, the propor¬ 
tion of CEMP stars formed from systems with initial secondary 
masses less than 0.7 Mq is 41% with model sets A and C, while 
it is 24% with model set B. 


In our model set D we inhibit the thermohaline mixing, con¬ 
sequently the material accreted by the secondary stars remains 
on the surface until it is mixed in by convection. As a result, 
the range of initial periods in which CEMP stars are formed in¬ 
creases (by roughly a factor of 25) because even a tiny amount 
of mass transferred from the AGB primary star on top of a dwarf 
or turnoff star is sufficient to enhance carbon. Consequently, the 
fractions Pc.to and Pc.dwa both increase, as does the overall 
CEMP/VMP ratio. 

Our knowledge of the initial properties of the halo stellar 
population is very uncertain and difficult to constrain. To deter¬ 
mine how robust our results are to changes in the age, the IMP, 
and the initial distributions of separations and mass ratios of 
our models, we simulate populations of VMP binary stars with 
a variety of model sets that adopt different assumptions. The 



Fig. 4. Initial period distribution computed with model sets A, B, 
and C (solid, dot-dashed and dashed line, respectively). The 0.2 
dex bins are equally spaced in log[o(Pi/days). 


CEMP/VMP ratio typically varies between 7% and 12%, with 
the exception of model Ilo 5 , while the proportion of different 
stellar types among CEMP stars is roughly constant. Below we 
briefly summarise the results obtained with the different model 
sets. 


Variations in the IMF (model sets /koi ond /los)- With the 
IMP proposed by |Kroupa ( 2001) l primary stars above one so¬ 
lar mass are favoured compared to the default IMP because 
the IMP slope is slightly shallower; therefore, with model set 
Ikoi there is a small increase in Fq. Compared to the default 
IMP, the IMP suggested by |Lucatello et al.| ( |20051 l produces 
more stars with M > 0.9 Mq (which can contribute to the for¬ 
mation of CEMP stars) than stars with M < 0.9 Mq (which 
do not form CEMP stars), therefore Ilos is the model that 
predicts the largest CEMP/VMP ratio. 

Variations in the initial mass-ratio distribution (model 
sets Qp4 and Q1). Model set Qp4, with <p{q) cx favours 
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MlogP 


Fig. 5. Left: Fq calculated with model sets that assume a power-law initial mass-ratio distribution, (p(q) oc with different values 
of Right: Fq of model sets that assume a log-normal initial-period distribution with different values of the mean, piogp, and fixed 
standard deviation, cr\ogp — 2.3. The filled circles show the value of Fq computed with our default model set A. 


systems with low mass ratios, q - Binary systems 

with low-mass secondary stars (M 2 _i < 0.5 Mq) do not con¬ 
tribute significantly to the CEMP fraction, either because the 
secondary stars need to accrete large amounts of mass, or 
they are not sufficiently luminous to be detected. As a con¬ 
sequence, Fc decreases. On the contrary, model set Ql, with 
fiq) oc q, favours equal mass ratios; therefore, binary sys¬ 
tems with relatively massive secondary stars are more com¬ 
mon, which need to accrete a small amount of mass to be¬ 
come sufficiently luminous CEMP stars. Consequently, Fq 
increases. In the left panel of Pig.|^we show the value of Fq 
computed with hve model sets that assume different y^ in a 
power-law distribution of the initial mass-ratio, fiq) oc g'A. 
The maximum CEMP/VMP ratio is found with = 1, i.e. 
with model set Ql. 

Variations in the initial-period distribution (model sets 5dm 
and S3). Because the initial-period distribution of halo bi¬ 
nary stars is poorly constrained, we test a variety of mod¬ 
els with log-normal period distributions in which we vary 
the mean, piogp, while we adopt a fixed standard deviation, 
crlog p - 2.3, i.e. the value determined by Duquennoy & 


|Mayor| fl991| l for a solar-neighbourhood population of bi- 
naries with G-dwarf primaries. The CEMP/VMP ratio com¬ 
puted with these models is shown in the right panel of Fig. 
as a function of piogp. The maximum is found with model 
set Sdmj that is, the Duquennoy & Mayor ( |1991| l distribu¬ 
tion with piogp = 4.8. With model set Sdm, the initial pe¬ 
riod distribution peaks at orbital periods between 5 x 10"^ 
and 10^ days, which is approximately the range where our 
WRLOF model of wind accretion is most efficient in form¬ 
ing CEMP stars. Because the average orbital period in the 
sample of 15 observed CEMP-i stars analysed in Paper I 
is approximately 1,500 days, we assume piogp = 3 in our 
model set S3. The CEMP/VMP ratio decreases compared 
to model Sdm because model S3 predicts more close binary 
systems that undergo a common-envelope phase and do not 
form CEMP stars. On the other hand, model set S3 predicts a 
larger CEMP/VMP ratio than our default model set A (with 
a flat log P-distribution) because it has fewer very wide sys¬ 
tems that do not interact. If a log-normal period distribution 
(as in model Sdm) is combined with a mass-ratio distribution 


that favours equal masses (as in model Ql), a CEMP/VMP 
fraction of approximately 15% is produced, i.e. almost the 
same fraction as assuming an IMF weighted towards rela¬ 
tively large masses, M > Q.9Mq. 

Variations in the minimum age of the synthetic population 
(model sets T8 and T12). A change in the minimum age of 
the stellar population in the halo by +2 Gyr has a small effect 
on Fq, as also found by Izzard et al. ( |2009 1 . This is mainly 
because of the difference in lifetime between the stars that 
accreted mass during the evolution and the stars that did not. 
As we discuss in more detail in Sect. |33] for a minimum 
age of the halo of fmin = 10 Gyr, stars of mass Q.9Mq and 
above have already become white dwarfs, whereas CEMP 
stars that (after accretion) have masses up to about 0.95Mq 
may still be visible. That is, only CEMP stars can be detected 
in the mass range [0.9,0.95]Mo for this age (Fig. 0. With 
fmin = 8 Gyr, CEMP stars up to about 1 .05Mq are still visible, 
whereas normally at this age stars more massive than Q.95Mq 
have already become white dwarfs, that is, the difference in 
the maximum mass of visible stars increases. As a result, 
the CEMP fraction increases. On the other hand, with = 
12 Gyr such a difference in mass decreases (the maximum 
masses are approximately O. 88 M 0 and 0.91Mo for normal 
stars and CEMP stars, respectively), and so does the CEMP 
fraction. 


3.4. Orbital periods of the CEMP stars. 

We now compare the orbital-period distribution predicted by 
our models with the observed sample of 15 binary CEMP stars 
with measured periods studied in Paper I. We exclude from this 
sample the stars CS 22964-161 A and B because they have 
most likely been polluted by a third star whose period is un¬ 
known ( [Thompson et al.|2008| l, LP 625-55 because it only has 
a lower limit on the orbital period, and CS 29497-034 which 
has [Fe/H] = -2.96. Our models produce most CEMP stars 
in orbits longer than a few thousand days, whereas the ob¬ 
served systems have orbital periods shorter than 4,600 days. 
However, in the SAGA sample, only 11 out of 94 CEMP stars 
with -2.8 < [Fe/H] < -1.8 have measured orbital periods, thus 
they may not be representative of the period distribution of the 
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Pf/days 



Fig. 6. Cumulative orbital-period (Pf) distributions for SAGA 
CEMP stars with -2.8 < [Fe/H] < -1.8 (shaded area). The 
lower limit of the observed distribution is determined assuming 
that all CEMP stars are binaries, whereas to determine the up¬ 
per limit we assume that only CEMP-i stars are in binaries. Our 
model sets A and Ilos are shown as solid and dotted lines, respec¬ 
tively. The 0.2 dex bins of the model distributions are equally 
spaced in logjQ(Pf/days). 


entire sample. To compare the observations with the results of 
our models, Fig.|^shows the observed and modelled cumulative 
period distributions. The implicit assumptions made in deriving 
the observed cumulative distribution are (/) that all CEMP stars 
are formed in binaries, most of which have unknown periods, 
and (ii) that the SAGA sample is complete for Porb ^ 4,600 days. 
The first assumption is probably correct for CEMP-i stars, which 
are mostly found in binary systems, but it is unlikely to be valid 
for CEMP-no stars (|Lucatello et al.||2006[ [Norris et al.||201^ 
Starkenburg et al.|2014[ Hansen et al. in prep.). For this reason 
we plot the cumulative period distribution as a range (shaded 
area in Fig. |^, in which the lower limit corresponds to the hy¬ 
pothesis that all 94 CEMP stars in the SAGA sample are binaries, 
while the upper limit corresponds to the hypothesis that only the 
71 CEMP-i stars in the sample are binaries. The solid and dotted 
lines in Fig. [^represent the cumulative period distributions pre¬ 
dicted with models A and Ilos, respectively, which bracket the 
distributions found with all the other models. 

Under the above assumptions, our models are consistent with 
the observations in the range between approximately 600 and 
5,000 days, and they seemingly predict an excess of CEMP stars 
with periods between a few tens and a few hundreds of days. 
On the other hand, the models do not reproduce the single ob¬ 
served CEMP staij^at Pf < 10 days. However, the assumption 
that the unknown periods of all other CEMP stars in our sample 
are longer than about 5,000 days may be incorrect; in reality at 
least some of these stars are likely to have shorter periods (e.g. 

. In addition, the comparison in Fig.|6 is 
based on the implicit assumption that our models reproduce the 
observed CEMP fraction. However, the CEMP/VMP ratio deter¬ 
mined for the SAGA sample is larger than predicted by a factor 
of 1.5 up to 3, depending on the model set (Table |^. Hence, 
given all the uncertainties, the results of this comparison are in¬ 
conclusive. 


Starkenburg et al.|2014 1 


^ Possible formation scenarios for this peculiar CEMP star with pe¬ 
riod 3.14 days, HE 0024-2523, are discussed, for example, byjLucatelloj 


et al. 


( 2003^ and in Paper I. 


Fig. 7. Mass distribution of CEMP stars (solid line) and C- 
normal stars (dot-dashed line) computed with model set A. Bins 
of width 0.025 dex are equally spaced in M»/Mo. 


3.5. Masses of CEMP stars 

Our models predict that approximately 90% of the CEMP stars 
have masses between 0.75 Mq and 0.95 Mq, with the peak of 
the distribution around 0.85 Mq. In Fig. we compare the mass 
distributions computed with model set A for CEMP and carbon- 
normal stars (solid and dot-dashed lines, respectively). The dis¬ 
tribution of CEMP stars is shifted towards higher mass by about 
0.05 Mq compared to the distribution of carbon-normal stars. 
This difference is a consequence of the mass-accretion process. 
A single star of mass 0.85 Mq ascends the giant branch after 
approximately 10.8 Gyr, and about 0.4 Gyr later it becomes 
a white dwarf. By contrast, if a star with an initial mass of 
0.75 Mq accretes 0.1 Mq from a binary companion and becomes 
a CEMP star, its life is longer than a normal 0.85 Mq star be¬ 
cause until mass accretion occurs the CEMP star evolves as a 
lower-mass star, i.e. more slowly. Currently, no mass measure¬ 
ments of CEMP stars are available for comparison with our mod¬ 
els, but future observations that aim to detect very metal-poor 
stars in eclipsing binary systems may allow the determination of 
whether the predicted difference in mass is correct. 


3.6. Abundance distributions 


The SAGA database includes the observed abundances of many 
elements. We compare the abundance distributions derived from 
the observations with those of our models. This comparison is 
qualitative because the observed sample is inhomogeneous and 
incomplete, but it can provide indications to constrain our mod¬ 
els. The magnitude limits adopted in the models are Vniin = 6 
and Umax = 16.5, consistent with the magnitude range of the very 
metal-poor stars in the observed sample. In the first part of this 
section we study the consequences of adopting different masses 
of the partial mixing zone (Sect. 3.6.11 while in the second part 


we compare the predictions of different model sets (Sect. 3.6.2 1 . 


3.6.1. The effect of Mpmz 

The mass of the partial mixing zone is a free parameter in our 
models. To understand its impact on the chemical properties of a 
population of CEMP stars we compare the element distributions 
computed with default model set A and four different values of 
MpMz- no PMZ, and Mpmz equal to 2 x 10“"^ Mq, 2 x 10“^ Mq 
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Table 4. Fractions of CEMP stars derived from the SAGA database and computed with our models. The adopted V-magnitudes 
range in the models is [6, 16.5]. CEMP stars have carbon abundance [C/Ee] > +0.7. 


model sets and data 

Fc (%) 

Cc,,.a(%) 

F C.TO (%) 

Cc.dwa (%) 

F gia.CEMP (%) 

F TO,CEMP ('%’) 

F dwa.CEMP (%) 

SAGA database 

27.98 ± 2.56 

22.56 ± 2.66 

38.95 + 7.55 

41.67+12.80 

62.86 + 8.55 

26.43 + 4.89 

10.71 + 2.91 

A 

9.47 ± 0.01 

10.22 ± 0.02 

10.79 + 0.04 

7.49 + 0.02 

64.25 + 0.12 

12.17 + 0.04 

23.58 + 0.06 

B 

9.95 ± 0.01 

10.85 + 0.02 

10.90 + 0.04 

7.83 + 0.02 

64.70 + 0.12 

11.70 + 0.04 

23.60 + 0.06 

D 

11.67 + 0.04 

10.57 + 0.05 

15.04 + 0.14 

12.64 + 0.08 

53.94 + 0.30 

13.77 + 0.13 

32.29 + 0.22 

Il05 

17.24 ± 0.05 

17.80 + 0.07 

18.26 + 0.16 

15.49 + 0.09 

64.04 + 0.29 

11.91+0.10 

24.05 + 0.15 



[hs/Fe] [Pb/Fe] 


Fig. 8. Panels (a)-(d): abundance distributions of sodium, light-i elements, heavy-i elements and lead, respectively. The 0.2 dex 
bins are equally spaced in [El/Ee]. The histograms represent CEMP and CEMP-s/r stars from our database (light and dark grey, 
respectively), with Poisson errors. Models with Mpmz = 0,2 x 10“"^ Mq, 2 x 10“^ Mq, 4 x IO^^Mq are shown as dotted, dashed, 
solid, dot-dashed lines, respectively. 


and 4 x 10“^ Mq. Pigs.[^-d show the abundance distributions of 
sodium, light-i elements, heavy-i elements]^ and lead. The his¬ 
tograms represent the distributions of the observed CEMP stars 
and the subsample of CEMP-s/r stars (light and dark grey, re¬ 
spectively). CEMP-s/r stars are plotted separately because their 
history of nucleosynthesis is probably different from CEMP stars 


with low abundances of europium and r-elements (e.g. Jonsell 

et al. 

20061|Sneden etal.|20081 

Lugaro et al.|2009l Bisterzo et al. 

2012 

lAbate et al.||2015a|b|l. Our models do not include the r- 


The abundances of light-i and heavy-J elements are de¬ 
fined as [Is/Fe] = ([Sr/Fe]+[Y/Fe]+[Zr/Fe])/3 and [hs/Fe] = 
([Ba/Fe] + [La/Fe]+[Ce/Fe])/3, respectively. If one of the elements 
is missing, we average the abundances of the two available elements. 


process, therefore in CEMP-i/r stars we do not expect to repro¬ 
duce the abundances of elements with a strong r-process com¬ 
ponent. However, the r-process is responsible for the production 
of many neutron-capture elements, including some amounts of 
elements traditionally associated with the s-process. In Paper II 
we show that CEMP-i/r stars have different abundance distri¬ 
butions than r-normal CEMP-i stars (i.e. with [Eu/Ee] < +1). 
In particular, the abundances of neutron-rich elements such as 
barium, lanthanum, cerium, and lead are typically strongly en¬ 
hanced and the [hs/ls] ratio is larger than is observed in r-normal 
CEMP-s stars and than predicted by our models, in line with 
the results of previous studies (e.g. |Bisterzo et al.|2012HLugaro 
let al.|2012D . 
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Fig. 9. Panels (a)-(d): distributions of [C/Fe], [Sr/Fe], [La/Fe], and [Eu/Fe], respectively. Model sets A, B, D, and Ilos are shown 
with solid, dot-dashed, dotted, and dashed lines, respectively. Light and dark grey histograms represent the observed CEMP and 
CEMP-s/r stars in the SAGA database. The 0.2 dex bins are equally spaced in [El/Ee]. 


In the models with a small PMZ (Mpmz = 0 and 2x10“"^ Mq), 
the majority of CEMP stars have normal, non-enhanced abun¬ 
dances of sodium (Pig. [^). This is because the abundance of 
sodium produced in AGB nucleosynthesis is sensitive to both the 
mass of the star and of its PMZ. With no PMZ, or only a small 
PMZ, only relatively massive AGB stars produce significant 
amounts of sodium. In the modelled population, most CEMP 
stars are formed in systems with a low-mass primary star, hence 
the distribution peaks around the initial abundance. On the other 
hand, with a relatively large PMZ (Mpmz = 2-4x 10“^ Mq), low- 
mass AGB stars also produce signihcant amounts of sodium; 
consequently the modelled distribution qualitatively reproduces 
the observations. 

The distributions of light-i and heavy-i elements (Pigs. 
and[^) do not vary significantly as a function of Mpmz, because 
low-mass AGB stars (below ^ 2 Mq) experience episodes of pro¬ 
ton ingestion from the convective envelope in the He-flash in¬ 
duced convection zone ( |Lugaro et al.|2012) l, and some amounts 
of s-elements are produced regardless of the PM^ Because 
low-mass primary stars dominate our synthetic binary popula¬ 
tion, the distributions of light-i and heavy-i elements in CEMP 


^ The nucleosynthetic signature of a proton-ingestion event is un¬ 
certain because one-dimensional stellar-evolution models do not ade¬ 
quately treat the physics of these inherently three-dimensional phenom¬ 
ena (see e.g. ; 


Stancliffe et al. 


2011 


Herwig et al. 2014 1 . 


stars are quite similar. By contrast, the abundance of lead is re¬ 
lated more strongly to the PMZ mass, as shown in Pig. |^, and 
the distributions obtained with small and large Mpmz are shifted 
by about 0.5 dex. Models with Mpmz > 2 x 10“^ Mq are in 
better agreement with the observed distribution; however, none 
of the models can reproduce the 11 observed CEMP stars with 
[Pb/Pe] > h- 2.8. Nine of these stars are CEMP-i/r stars, and the 
discrepancy with the models may indicate that lead in these ex¬ 
tremely enriched stars was not produced purely by the s-process. 
We note that the observed abundances are determined assuming 
local thermodynamic equilibrium. It has been shown that depar¬ 
tures from this approximation can greatly impact the abundances 
of giants and metal-poor stars (e.g. [Bergemann et al.||2012| l. In 
particular, positive corrections of 0.3-0.5 dex are to be expected 
for the lead abundances of stars in the metallicity range of our 
observed sample ( Mashonkina et al.|2012| l. If we assume an av¬ 
erage correction of 0.4 dex, 21 CEMP-s stars (of which 16 are 
CEMP-i/r) would have [Pb/Pe] > 2.8 and cannot be reproduced 
by any of our models. 


The results of the comparison between models with different 
partial mixing zones, as well as the results of Papers I and II, 
suggest that to reproduce the chemical properties of the observed 
sample our models require a relatively large Mpmz- Therefore, in 
the following we adopt Mpmz = 2 x 10 “^ Mq in our model sets. 
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3.6.2. Abundance distributions with different model sets 


Our modelled abundance distributions depend mostly on the dis¬ 
tributions of the primary and secondary stellar mass of our syn¬ 
thetic CEMP stars; the former determines the amount of each el¬ 
ement produced by AGB nucleosynthesis, while the latter essen¬ 
tially determines how much mass is accreted and consequently 
the abundance enhancement. We compare the results of our de¬ 
fault model set A with three model sets: model set B, in which 
the M 2 _i distribution of CEMP stars is more peaked towards rel¬ 
atively massive secondary stars (M 2 ,i > 0.7 Mq), model set Ilos, 
in which a significant proportion of CEMP stars are formed in 
binary systems with a primary star initially above 1.2 Mq, and 
model set D, in which thermohaline mixing is inhibited and 
therefore the abundances of the accreted material are not di¬ 
luted in turnoff and dwarf CEMP stars. The CEMP fractions 
observed in the SAGA database and computed with these four 
model sets are shown in Table In model sets A, B, and 1 lo5 , 
the CEMP/VMP ratios are higher than the corresponding sets in 
Table [^because for the comparison with SAGA stars we adopt 
lower magnitude limits in our simulations, thus the fainter dwarf 
stars are selected against. These fainter stars have a larger pro¬ 
portion of carbon-normal stars, so the CEMP/VMP ratios in¬ 
crease. On the other hand, in model set D, among the fainter 
dwarf stars the CEMP stars are more frequent because of the 
lack of dilution, thus Fq decreases. 


Eigures [^- ^ show the abundance distributions of carbon, 
strontium, lanthanum, and europium. We find similar trends 
for all elements. Compared to the default model A, model set 
Ilo 5 favours more massive primary stars, which typically pro¬ 
duce larger abundances of carbon and i-elements for Mpmz = 
2x10^^ Mq. Also, in model set 1 lo5 more CEMP stars are formed 
with M 2 ,i < 0.6 Mq, which accrete large amounts of material 
from the primary star. Consequently, the abundance distributions 
of model set Ipos are weighted towards large enhancements. On 
the contrary, model set B favours relatively massive secondary 
stars, M 2 ,i > 0.7 Mq, which accrete small amounts of material 
that is more strongly diluted. Consequently, model set B predicts 
abundance distributions that peak at lower [El/Ee] than our de¬ 
fault model A. With model set D, dwarfs form a clearly distinct 
group of stars. Because the accreted material is not efficiently 
mixed throughout the entire star, these stars have a chemical 
composition very similar to that of the polluting stars, and hence 
the largest enhancements. 

The distributions of lanthanum and europium (respectively, 
a heavy-i and an r-element) in Pigs. and show that the 
abundances observed in CEMP-i/r stars are outside the range 
that can be produced by our models. If we exclude the CEMP- 
s/r stars from the sample, our models qualitatively reproduce the 
observed distributions. This discrepancy is observed also in the 
distributions of [hs/Pe] and [Pb/Pe] (Pigs.|^ and[^). By con¬ 
trast, the distributions of carbon, sodium, strontium, and light-i 
elements do not differ significantly between CEMP-s/r stars and 
r-normal CEMP-s stars. These results suggest that there is a re¬ 
lation between the strongly enhanced abundances of the heavy- 
s elements and the r-elements in the population of CEMP-i/r 
stars, in line with our conclusions in Paper 11. 

The origin of the r-enhancement in CEMP-s/r stars is un¬ 
certain. If we start our simulations with an initially enhanced 
abundance of r-elements to mimic the effect of a primordial en¬ 
richment independent of the abundances of i-elements, as sug¬ 
gested e.g. by Bisterzo et al. ( 2012[ ), our models fail to repro¬ 
duce the stars with the lowest [Eu/Fe] without improving the 
results obtained for CEMP-s/r stars. Pig. [T0| shows the observed 



[Eu/Fe] 


Fig. 10. Distribution of [Ba/Pe] vs. [Eu/Pe]. The red distribution 
represents our default model A. The grey distribution is com¬ 
puted with an initial enhancement of [r/Pe]ini = 1. The dotted 
lines indicate the thresholds [Ba/Pe] = -1-0.5 and [Eu/Pe] = 
1 that define CEMP-s and CEMP-s/r stars, respectively. The 
crosses indicate the observed CEMP-s stars with [Eu/Pe] < -i-l. 
CEMP-s/r are shown as blue circles. The -H symbol in the top- 
left corner shows the average observed uncertainty. 


abundances of europium and barium in the CEMP stars of our 
sample, together with the synthetic abundances computed with 
our default model set A (red distribution), assuming an initial 
enhancement of [r/Pe]ini = 1 (grey distribution). The observa¬ 
tions exhibit a clear linear correlation between the abundances 
of barium and europium. Our model set A qualitatively repro¬ 
duces this correlation for r-normal CEMP-i stars, that is, to the 
left of the vertical dotted line at [Eu/Pe] = -i-l. As discussed by 
Lugaro et al. ( 2012| ), an initial enhancement of [r/Pe]ini < 0.4 
dex is essentially washed out by the europium produced by the s- 
process in the AGB nucleosynthesis. If we adopt [r/Pe]ini = -i-l, 
europium increases while barium does not change, and the ob¬ 
served correlation is not reproduced. If we adopt [r/Pe]ini = +2 
to reproduce the most r-rich CEMP-s/r stars, then the abun¬ 
dances of the r-elements in our synthetic CEMP stars all end 
up in one bin around this value. 

Pigure [TT] shows the distribution of [Ba/Pe] vs. [C/Pe] com¬ 
puted with model set A. An apparent correlation exists between 
carbon and barium in the observed CEMP-s stars, regardless of 
the r-enhancement, although there is considerable scatter in the 
data. In the binary scenario for the formation of CEMP stars this 
correlation is a direct consequence of nucleosynthesis in AGB 
stars. The amounts of carbon and barium produced by an AGB 
star are correlated, and depend only on the stellar mass (for a 
fixed PMZ). The observed correlation is qualitatively reproduced 
by our model, although a fraction of the CEMP stars have a bar¬ 
ium abundance systematically larger than the model predictions. 
This discrepancy is not removed even with higher masses of the 
PMZ, and suggests that, for a given amount of carbon, our AGB 
models should produce a larger amount of barium. 

The observed correlation between carbon and barium breaks 
down for CEMP-no stars, whose carbon abundances are spread 
over two orders of magnitude, and are not well reproduced by 
our models. This discrepancy may arise within the mass-transfer 
scenario because a large abundance of carbon, combined with 
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Fig. 11. As Fig. 10 for [Ba/Fe] vs. [C/Fe]. The vertical line rep¬ 


resents the threshold [C/Fe] = h- 0.7. In the bottom-right corner 
the average observed uncertainty is shown as a -H symbol. 


a small abundance of barium, can be obtained in models of 
massive AGB stars (M > 2.8Mo) with no PMZ, and massive 
stars are relatively rare in our population. However, about half 
of the CEMP-no stars in our sample have [Ba/Fe] < 0, and 
with our model it is not possible to reproduce their abundances. 
The formation scenario of CEMP-no stars is currently uncertain. 
Masseron et al. (2010|l show that there is a clear continuity from 


CEMP-i to CEMP-no in the abundance trends of several ele¬ 
ments, including C-hN, O, Mg, '^C/^^C and C/N, pointing to¬ 
wards the scenario of mass transfer from an AGB star. On the 
other hand, the study of Starkenburg et al. (2014 1 indicates that 
CEMP-no stars are not found in binary systems more often than 
carbon-normal stars or, alternatively, they belong to very wide 
binaries. Although the hypothesis of wide separations is consis¬ 
tent with our results, the evidence that with our models most 
of the CEMP-no stars are not reproduced in the carbon-barium 
space may be an indication that their observed carbon enhance¬ 
ment is not related with AGB nucleosynthesis but rather is the 
result of contributions from first-generation stars. 


4. Discussion 

In this paper we have simulated populations of binary stars at 
metallicity Z = 10“"^ ([Ee/H ^ -2.3) to investigate the mass- 
transfer scenario for the formation of CEMP stars. A compari¬ 
son with similar studies (e.g. |Izzard et al.|2009| and |Abate~et al.| 
2013| l shows that our default model set A predicts a CEMP/VMP 
ratio more than three times larger than previous models adopt¬ 
ing default initial assumptions. In the default model adopted in 
[Abate et al.| ( |2013| l, AGB stars of mass below 1.2 Mq do not 
undergo third dredge-up, whereas in our models (based on the 
detailed calculations of |Karakas]|20 1 Oj ) the minimum mass for 
third dredge-up in AGB stars is 0.9 Mq. Therefore, the range of 
primary masses that contribute to the formation of CEMP stars 
is larger and, as a consequence, the CEMP/VMP ratio predicted 


in our models increases. In our simulations we adopt a selection 
criterion based on the luminosity of the stars, instead of a simple 
cut-off in surface gravity as in Abate et al. ( |2013 i. As a result, 
the proportion of dwarf stars in the CEMP population increases 
to approximately 50%, but the total CEMP/VMP ratio remains 
essentially the same. 

If we adopt a solar-neighbourhood IMP, our CEMP/VMP 
ratio, Fq, is approximately 8 % for [C/Ee] > -i-l, consistent 
with the observed value of the SEGUE sample in the range 
-2.5 < [Pe/H] < -2.0. Therefore, a change of the IMP towards 
relatively large primary masses (as suggested, e.g. by Lucatello 
et al. [2005 and Suda et al. [201 3|l is not necessary to repro¬ 


duce the observed CEMP/VMP ratio. However, the proportion 
of dwarf CEMP stars predicted in our models is much larger 
than observed in the SEGUE sample at -2.5 < [Pe/H] < -2.0. 
Also, the total number of dwarf stars observed at this metallicity 
in the SEGUE sample is a factor of four to five lower than the 
number of turnoff and giant stars. These differences may suggest 
that our simulations include stars that in reality are too faint to 
be detected, or that the observed surface gravities are possibly 
underestimated by a few tenths of a dex. 

The CEMP/VMP ratio observed in the SEGUE sample in¬ 
creases by almost a factor of two if the minimum carbon abun¬ 
dance in the definition of CEMP stars is reduced from [C/Pe] = 
-H 1.0 to [C/Pe] = h- 0.7. By contrast, this causes only a small dif¬ 
ference in our models. In Paper II we showed that most carbon- 
normal VMP stars from the SAGA database exhibit carbon abun¬ 
dances within approximately 0.4 dex from the average value of 
[C/Pe] = h- 0.3. Assuming that carbon-normal VMP stars in 
the SEGUE sample have the same distribution of carbon abun¬ 
dances, and considering that an observational error of 0.3 dex 
is associated with carbon abundances in SEGUE stars, one ex¬ 
planation of the large increase in the observed CEMP fraction 
between [C/Pe] = -1-0.7 and [C/Pe] = H-1.0 is that we are ob¬ 
serving the tail of the carbon-abundance distribution of carbon- 
normal very metal-poor stars. 

The CEMP/VMP ratio calculated among stars in the SAGA 
database is very different from the SEGUE sample. In the 
range -2.8 < [Pe/H] < -1.8 we find Fq - 28%, adopting 
[C/Pe]njin = - 1 - 0 . 7 , almost a factor of three larger than the value 
calculated for SEGUE stars. A discrepancy of about a factor of 
three is also found with most of our model sets, which predict 
CEMP/VMP ratios around 10%. Another remarkable difference 
with the SEGUE sample is that the CEMP fraction from the 
SAGA sample increases towards stars of higher logjg g. A possi¬ 
ble explanation of these discrepancies is that the SAGA sample 
is a compilation of all very metal-poor stars currently available 
in the literature. Hence, it is an inhomogeneous and incomplete 
sample, and is possibly biased towards chemically peculiar stars. 
However, the large differences between the results obtained with 
the SEGUE and SAGA samples point out the general problem 
that the measured fractions depend significantly on the observed 
sample considered. Consequently, the results of the comparisons 
with population synthesis models are as yet inconclusive. 

In our simulations we assume a binary fraction of unity in 
the range of orbital separations between 5Rq and 5 x 10^ Rq. 


Eor comparison, Kouwenhoven et al. (2007 1 find that, with 3cr 
confidence, the binary fraction in the young stellar association 
Scorpius OB2 is at least 0.7, and probably close to unity for 
stars of spectral type A and B in the same range of separa¬ 
tions. However, the binary fraction of the stellar population of 
the Galactic halo is poorly constrained ([Carney et al. 2001 2005 [ 


Rastegaev|2010[ Gao et al.|2014 Aoki et al.|2015[ ). Recent re¬ 

sults suggest that among stars with [Ee/H] < -1.1 the binary 
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fraction is at least 50% and possibly increasing with decreasing 
metallicity ( Gao et al.|2014[ Yuan et al.|2015 i 


All our model sets predict that CEMP stars are formed in ini¬ 
tially wide binary systems with orbital periods typically longer 
than 1,000 days up to approximately 10^ days. About 80% of 
the simulated CEMP stars have periods longer than 5,000 days, 
which is approximately the longest period currently observed. 
The median of the observed period distribution is approximately 
500 days, consistent with the average orbital period estimated 
by IStarkenburg et al.| ( |2014l for a different sample of CEMP-s 
stars, whereas for the synthetic distributions it is about 20,000 
days. This indicates that, if the distribution determined from the 
CEMP stars with measured orbital periods is representative of 
the entire CEMP population, our simulations signihcantly un¬ 
derestimate the observed CEMP/VMP ratio. A similar discrep¬ 
ancy between the observed periods and model predictions also 
exists at higher metallicity for the barium stars, which are likely 
formed by the same wind-accretion process as CEMP stars (e.g. 
[Boffin & Jorissen|1988[|Jorissen et al.|1998| l. Eor barium stars, 
the observed sample of orbital periods is complete, and only 
a few systems are known with periods longer than lO'* days. 
The discrepancy between the observed and modelled periods 
suggests that in our models binary systems need to efficiently 
lose angular momentum, and transfer mass more efficiently at 
short separations, as also discussed by jlzzard et al.| (20 T0] l and in 
Papers I and II. 

The comparison between the observed and synthetic abun¬ 
dance distributions provides constraints on our nucleosynthe¬ 
sis model and on the properties of the progenitor systems. The 
chemical composition of a modelled CEMP star depends essen¬ 
tially on the mass of the primary star, which determines the ele¬ 
ments produced by AGB nucleosynthesis, and on the initial mass 
of the secondary star, which determines the amount of material 
that needs to be accreted so that after 10 Gyr the star is luminous 
enough to be visible and enriched in carbon. Most model sets 
predict similar distributions of Mi_i and M 2 ,i, consequently the 
abundance distributions are similar. 

The modelled abundance distributions of most elements do 
not show a strong dependence on the mass of the partial mix¬ 
ing zone, MpMz- The production of s-elements is sensitive to 
MpMz, especially for stars with > 1.5 Mq, while most 
CEMP stars in our simulations are formed from binary systems 
with Mi i < 1.2 Mq. By contrast, in Papers I and II we hnd 
that primary masses above 1.4 Mq are necessary to reproduce 
the detailed chemical composition of most observed CEMP-s 
stars of our sample. Systems with primary mass between 1.5 Mq 
and 2.5 Mq are required to reproduce the element-to-element 
ratios observed in CEMP-i stars that exhibit abundant heavy- 
s elements and large [hs/ls] and [Pb/hs] values. These chemi¬ 
cal properties are observed mostly in CEMP-s/r stars, and also 
in some r-normal CEMP-s stars. The discrepancy between the 
results of the population-synthesis simulations and the detailed 
analysis of Paper I and II may hint that nucleosynthesis models 
of low-mass AGB stars (M <1.2 Mq) should produce, in some 
circumstances, higher abundances of heavy-i elements and lead. 

We note that the observed abundances are normally deter¬ 
mined under the assumption of local thermodynamic equilib¬ 
rium, which is not always fuffilled in the layers of the stellar 


atmosphere where spectral lines are formed (e.g. Mihalas & 
Athay |1973|l. The departure from local thermodynamic equi- 


stars should be applied to the observed lead abundances. This 
correction increases for increasing temperature, and for decreas¬ 
ing metallicity and surface gravity. Positive corrections are typi¬ 
cally required also for barium and europium ( [Mashonkina et ah] 
2012t Niemczura et al.||2014|l. Although the magnitudes of the 


corrections are of the order of the observational uncertainties 
(< 0.2 dex), they suggest that our models should produce higher 
abundances of heavy elements for a given amount of carbon. 

Our default model set A predicts abundance distributions of 
light elements (e.g. carbon and sodium) and light-s elements 
that are qualitatively consistent with the observed distributions 
in all CEMP stars. In contrast, the abundance distributions of 
heavy-i elements, lead, and europium are reproduced only in r- 
normal CEMP stars, whereas in CEMP-s/r stars the abundances 
of these elements are outside the range of our model predictions. 
Similarly, a clear correlation between the abundances of barium 
and europium is observed in all CEMP-i stars. Our model re¬ 
produces such a correlation in CEMP stars with [Eu/Ee] < -i-l, 
while in CEMP-i/r stars the observed enhancements of barium 
and europium are too large. Consistent results are found in the 
detailed analysis of CEMP-s stars performed in Papers I and 
II. In Paper I we hnd that to reproduce the large abundances 
of heavy-i elements observed in CEMP-i/r stars, our models 
typically overestimate the abundances of carbon, sodium, mag¬ 
nesium, and light-s elements. Paper II shows that in r-normal 
CEMP-s stars the abundances of most elements are reproduced 
on average within the observed uncertainty, whereas in CEMP- 
i/r stars the abundances of both heavy-i and r-elements are sys¬ 
tematically underestimated. 

Although enhanced abundances of r-elements are also ob¬ 
served in carbon-normal stars with no evidence of duplicity (e.g. 
[Barklem et al.|[2005] [Roederer et al.|[2014a|b[ Hansen et al. in 
prep.), our results suggest that the s- and r-processes responsi¬ 
ble for the abundances observed in CEMP-i/r stars may have 
occurred in the same astrophysical site. Consequently, these re¬ 
sults indicate that under some conditions AGB stars may be able 
to reach the large densities necessary to activate the r-process, 
in contrast to what is currently found in the AGB models. 
Detailed models of AGB stars at extremely low or zer o metallic 


ity show that large neutron densities (lO'^ - 10'®'’’^ ^ 


“cm 


Campbell 


& Lattanzio|2008| Herwig et al.|201 11 may be reached as a re- 


lihrium in real stars introduces a bias in the abundance deter¬ 
mination of most elements, including carbon, iron, strontium, 
barium, and lead. Eor example, [Mashonkina et ah (20121 show 
that a positive correction of about 0.3-0.5 dex for metal-poor 


suit of proton ingestion in the helium-flash region. [Lugaro et ah 
( 2QQ9\ propose a speculative scenario in which extremely low- 
metallicity (Z < 10“^) AGB stars may be able to produce both s- 
and r-elements. Eurther evolutionary and nucleosynthetic calcu¬ 
lations are needed to test quantitatively if this scenario may also 
work at metallicity Z ^ 10 and if the predicted abundances re¬ 
produce the chemical compositions of observed CEMP-i/r stars. 

5. Conclusions 

In conclusion, in our synthetic population of very metal-poor bi¬ 
nary stars we hnd a CEMP fraction that varies between 7% and 
17%, depending on the initial assumptions. Our default model 
set predicts a CEMP/VMP fraction of 8.5% for [C/Ee] > -1-0.7. 
This fraction is more than three times larger than the results ob¬ 
tained by Izzard et al. ( 2009| l and Abate et al. ( 2013[ ) with the 
same population-synthesis code and default initial assumptions. 
This difference results from the updates that we included in the 
model of AGB nucleosynthesis. With the default models adopted 
in the previous work, AGB stars below 1.2 Mq did not undergo 
third dredge-up, whereas in our updated models AGB stars ex¬ 
perience third dredge-up down to 0.9 Mq. Hence, the range of 
primary masses that can form a CEMP star is increased. 
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The observed fraction of CEMP stars depends significantly 
on the sample that is taken into account. In the most recent study 
of the SDSS/SEGUE stellar survey a CEMP fraction of 11.5% 
is found for stars with -2.5 < [Ee/H] < -2. Considering the 
large uncertainties associated at low metallicity with the ini¬ 
tial distributions of masses and separations in binary systems, 
the CEMP/VMP ratio predicted with our models is consistent 
with the observed value, although the models predict many more 
dwarf CEMP stars than are observed. 

There are currently few measured orbital periods for CEMP 
stars, and therefore it is difficult to provide strong constraints on 
our model of angular-momentum loss and wind-accretion effi¬ 
ciency. Euture observations of the orbital periods of additional 
CEMP stars will allow us to gain insight into the mass-transfer 
process at low metallicity. At present, the orbital periods of most 
of our synthetic stars are on average ten times longer than the 
typical observed periods. This indicates that our models should 
produce more CEMP stars in binary systems below a few times 
10^ days, or, alternatively, that most CEMP stars have periods 
longer than 10'^ days. 

The population of synthetic CEMP stars qualitatively repro¬ 
duces the abundance distributions of carbon, sodium, and light-s 
elements observed in all CEMP stars. On the other hand, the ob¬ 
served distributions of heavy-s elements, lead, and europium are 
reproduced only in r-normal CEMP-s stars, whereas in CEMP- 
sjr stars the abundances of these elements are underestimated. 
The correlation between the abundances of barium and europium 
observed in all CEMP stars indicate that the enhancements in s- 
and r-elements are not independent. This suggests that both the 
s- and r-processes may be active in the same astrophysical site. 

Acknowledgements. We thank the referee for her/his helpful and supportive 
comments on our paper. CA is grateful to P. van Oirschot for many use¬ 
ful discussions on the stellar distribution in the Galactic halo and to the 
Netherlands Organisation for Scientific Research (NWO) for funding support 
under grant 614.000.901. RJS is the recipient of a Sofja Kovalevskaja Award 
from the Alexander von Humboldt Foundation. RGI is grateful to the Alexander 
von Humboldt Foundation and the Science and Technology Facilities Council 
(STFC) for funding support. AIK is supported through an Australian Research 
Council Future Fellowship (FTl 10100475). TCB acknowledges partial sup¬ 
port for this work from grants PHY 08-22648; Physics Frontier Center/Joint 
Institute for Nuclear Astrophysics (JINA), and PHY 14-30152; Physics Frontier 
Center/JINA Center for the Evolution of the Elements (JINA-CEE), awarded by 
the US National Science Foundation. YS.L. acknowledges partial support by the 
2014 Research Fund of the Chungnam National University. 


References 

Abate, C., Pols, O. R., Izzard, R. G., & Karakas, A. I. 2015b, A&A, in press 
Abate, C., Pols, O. R., Izzard, R. G., Mohamed, S. S., & de Mink, S. E. 2013, 
A&A, 552, A26 

Abate, C., Pols, O. R., Karakas, A. L, & Izzard, R. G. 2015a, A&A, 576, Al 18 
Aoki, W., Beers, T. C., Christlieb, N., et al. 2007, Apl, 655, 492 
Aoki, W., Beers, T. C., Lee, Y. S., et al. 2013, AJ, 145, 13 
Aoki, W., Norris, I. E., Ryan, S. G., Beers, T. C., & Ando, H. 2002, ApJ, 567, 
1166 

Aoki, W., Suda, T., Beers, T. C., & Honda, S. 2015, AJ, 149, 39 
Asplund, M., Grevesse, N., Sauval, A. I., & Scott, P. 2009, ARA&A, 47, 481 
Barklem, P. S., Christlieb, N., Beers, T. C., et al. 2005, A&A, 439, 129 
Beers, T. C. & Christlieb, N. 2005, ARA&A, 43, 531 
Beers, T. C., Preston, G. W., & Shectman, S. A. 1985, AJ, 90, 2089 
Beers, T. C., Preston, G. W., & Shectman, S. A. 1992, AJ, 103, 1987 
Bergemann, M., Lind, K., Collet, R., Magic, Z., & Asplund, M. 2012, MNRAS, 
427, 27 

Bessell, M. S., Castelli, E, & Plez, B. 1998, A&A, 333, 231 
Bisterzo, S., Gallino, R., Straniero, O., Cristallo, S., & Kappeler, F. 2012, 
MNRAS, 422, 849 

Boffin, H. M. J. & Jorissen, A. 1988, A&A, 205, 155 

Campbell, S. W. & Lattanzio, J. C. 2008, A&A, 490, 769 

Carney, B. W., Aguilar, L. A., Latham, D. W., & Laird, J. B. 2005, AJ, 129, 1886 


Carney, B. W., Latham, D. W., Laird, J. B., Grant, C. E., & Morse, J. A. 2001, 
AJ, 122, 3419 

Carollo, D., Beers, T. C., Bovy, J., et al. 2012, ApJ, 744, 195 

Christlieb, N., Green, P. J., Wisotzki, L., & Reimers, D. 2001, A&A, 375, 366 

Christlieb, N., Schorck, T., Erebel, A., et al. 2008, A&A, 484, 721 

Cohen, J. G., Shectman, S., Thompson, L, et al. 2005, ApJ, 633, L109 

Duquennoy, A. & Mayor, M. 1991, A&A, 248, 485 

Eggleton, P. P. 1971, MNRAS, 151, 351 

Erebel, A., Christlieb, N., Norris, J. E., et al. 2006, ApJ, 652, 1585 
Gao, S., Liu, C., Zhang, X., et al. 2014, ApJ, 788, L37 
Gunn, J. E., Siegmund, W. A., Mannery, E. J., et al. 2006, AJ, 131, 2332 
Helmi, A. 2008, A&A Rev., 15, 145 

Herwig, F., Pignatari, M., Woodward, P. R., et al. 2011, ApJ, 727, 89 
Herwig, F, Woodward, P. R., Lin, P.-H., Knox, M., & Fryer, C. 2014, ApJ, 792, 
L3 

Izzard, R. G., Dermine, T., & Church, R. P. 2010, A&A, 523, A10+ 

Izzard, R. G., Dray, L. M., Karakas, A. L, Lugaro, M., & Tout, C. A. 2006, A&A, 
460, 565 

Izzard, R. G., Glebbeek, E., Stancliffe, R. J., & Pols, O. R. 2009, A&A, 508, 
1359 

Izzard, R. G., Tout, C. A., Karakas, A. L, & Pols, O. R. 2004, MNRAS, 350, 407 
Jonsell, K., Barklem, P. S., Gustafsson, B., et al. 2006, A&A, 451, 651 
Jorissen, A., Van Eck, S., Mayor, M., & Udry, S. 1998, A&A, 332, 877 
Juric, M., Ivezic, Z., Brooks, A., et al. 2008, ApJ, 673, 864 
Karakas, A. & Lattanzio, J. C. 2007, PASA, 24, 103 
Karakas, A. 1. 2010, MNRAS, 403, 1413 

Karakas, A. L, Lattanzio, J. C., & Pols, 0. R. 2002, PASA, 19, 515 
Kobayashi, C., Karakas, A. L, & Umeda, H. 2011, MNRAS, 414, 3231 
Komiya, Y., Suda, T., Minaguchi, H., et al. 2007, ApJ, 658, 367 
Kouwenhoven, M. B. N., Brown, A. G. A., Portegies Zwart, S. F, & Kaper, L. 

2007, A&A, 474, 77 
Kroupa, P. 2001, MNRAS, 322, 231 

Kroupa, P, Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545 
Lee, Y. S., Beers, T. C., Masseron, T, et al. 2013, AJ, 146, 132 
Lee, Y. S., Suda, T., Beers, T. C., & Stancliffe, R. J. 2014, ApJ, 788, 131 
Lucatello, S., Beers, T. C., Christlieb, N., et al. 2006, ApJ, 652, L37 
Lucatello, S., Gratton, R., Cohen, J. G., et al. 2003, AJ, 125, 875 
Lucatello, S., Gratton, R. G., Beers, T. C., & Carretta, E. 2005, ApJ, 625, 833 
Lugaro, M., Campbell, S. W., & de Mink, S. E. 2009, PASA, 26, 322 
Lugaro, M., Karakas, A. L, Stancliffe, R. J., & Rijs, C. 2012, ApJ, 47, 1998 
Mackey, J., Bromm, V., & Hernquist, L. 2003, ApJ, 586, 1 
Maeder, A., Przybilla, N., Nieva, M.-F, et al. 2014, A&A, 565, A39 
Marsteller, B., Beers, T. C., Rossi, S., et al. 2005, Nuclear Physics A, 758, 312 
Mashonkina, L., Ryabtsev, A., & Frebel, A. 2012, A&A, 540, A98 
Masseron, T., Johnson, J. A., Plez, B., et al. 2010, A&A, 509, A93 
Meynet, G., Ekstrdm, S., & Maeder, A. 2006, A&A, 447, 623 
Meynet, G., Hirschi, R., Ekstrom, S., et al. 2010, A&A, 521, A30 
Mihalas, D. & Athay, R. G. 1973, ARA&A, 11, 187 

Moni Bidin, C., Carraro, G., Mendez, R. A., & Smith, R. 2012, ApJ, 751, 30 
Niemczura, E., Smalley, B., & Pych, W. 2014, Determination of Atmospheric 
Parameters of B-, A-, F- and G-Type Stars: Lectures from the School 
of Spectroscopic Data Analyses, GeoPlanet: Earth and Planetary Sciences 
(Springer) 

Nomoto, K., Kobayashi, C., & Tominaga, N. 2013, ARA&A, 51, 457 
Norris, J. E., Yong, D., Bessell, M. S., et al. 2013, ApJ, 762, 28 
Placco, V. M., Frebel, A., Beers, T. C., & Stancliffe, R. J. 2014, ApJ, 797, 21 
Pols, 0. R., Izzard, R. G., Stancliffe, R. J., & Glebbeek, E. 2012, A&A, 547, A76 
Pols, O. R., Schroder, K.-R, Hurley, J. R., Tout, C. A., & Eggleton, P. P. 1998, 
MNRAS, 298, 525 

Pols, O. R., Tout, C. A., Eggleton, P. R, & Han, Z. 1995, MNRAS, 274, 964 
Rastegaev, D. A. 2010, AJ, 140, 2013 

Reimers, D. 1975, Memoires of the Societe Royale des Sciences de Liege, 8, 369 
Richard, O., Michaud, G., & Richer, J. 2002a, ApJ, 580, 1100 
Richard, O., Michaud, G., Richer, J., et al. 2002b, ApJ, 568, 979 
Roederer, 1. U., Cowan, J. J., Preston, G. W., et al. 2014a, MNRAS, 445, 2970 
Roederer, 1. U., Preston, G. W., Thompson, 1. B., et al. 2014b, AJ, 147, 136 
Sandage, A. 1972, ApJ, 178, 25 

Sneden, C., Cowan, J. J., & Gallino, R. 2008, ARA&A, 46, 241 
Spada, F., Demarque, R, Kim, Y.-C., & Sills, A. 2013, ApJ, 776, 87 
Stancliffe, R. J. 2009, MNRAS, 394, 1051 

Stancliffe, R. J., Dearborn, D. S. P, Lattanzio, J. C., Heap, S. A., & Campbell, 
S.W. 2011, ApJ, 742, 121 

Stancliffe, R. J. & Eldridge, J. J. 2009, MNRAS, 396, 1699 
Stancliffe, R. J. & Glebbeek, E. 2008, MNRAS, 389, 1828 
Stancliffe, R. J., Glebbeek, E., Izzard, R. G., & Pols, O. R. 2007, A&A, 464, L57 
Starkenburg, E., Shetrone, M. D., McConnachie, A. W, & Venn, K. A. 2014, 
MNRAS, 441, 1217 

Suda, T., Katsuta, Y, Yamada, S., et al. 2008, PASJ, 60, 1159 


16 



C. Abate et al.: Modelling the observed properties of CEMP stars using binary population synthesis 


Suda, T., Komiya, Y., Yamada, S., et al. 2013, MNRAS, 432, L46 

Suda, T., Yamada, S., Katsuta, Y, et al. 2011, MNRAS, 412, 843 

Thompson, I. B., Ivans, 1.1., Bisterzo, S., et al. 2008, ApJ, 677, 556 

Toonen, S. & Nelemans, G. 2013, A&A, 557, A87 

Torres, G. 2010, AJ, 140, 1158 

Umeda, H. & Nomoto, K. 2003, Nature, 422, 871 

Umeda, H. & Nomoto, K. 2005, ApJ, 619, 427 

van Oirschot, R, Nelemans, G., Helmi, A., et al. 2014, A&A, 298, 431 

Vassiliadis, E. & Wood, P. R. 1993, ApJ, 413, 641 

Windhorst, R. A., Burstein, D., Mathis, D. F., et al. 1991, ApJ, 380, 362 

Yanny, B., Rockosi, C., Newberg, H. J., et al. 2009, AJ, 137, 4377 

Yong, D., Norris, J. E., Bessell, M. S., et al. 2013, ApJ, 762, 27 

York, D. G., Adelman, J., Anderson, Jr., J. E., et al. 2000, AJ, 120, 1579 

Yuan, H., Liu, X., Xiang, M., et al. 2015, ApJ, 799, 135 


17 



